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Abstract 

The Slater orbitals are the natural basis functions in quantum molecular calculations. Three- 
center repulsion Coulomb-exchange integrals over Slater orbitals are evaluated analytically 
with arbitrary orbital exponents, first for linear conformation of the atomic centers. These 
integrals have been expressed as a linear combination of three-center one-electron overlap 
integrals, and those have been calculated using auxiliary functions in terms of one-electron 
auxiliary integrals. Only one infinite expansion has been introduced. The resulting integral 
converges to 20 decimal digits using about 25-30 terms. Hybrid-exchange three-center repul- 
sion integrals will be investigated next using this method, as well as triangular conformation 
of the centers. 

keywords: Slater orbitals; three-center integrals; two-electron repulsion integrals 



f Permanent address: Department of Physics, University of Rize, Turkey. 
This work has been supported by DFG and TUBITAK. 



1. Introduction 



Slater-type orbitals (STO) pQ are the natural basis functions in quantum molecular calculations. 
STOs represent well the electron density near the nucleus (cusp) and far from the nucleus (correct 
asymptotic decay). Therefore the STOs resemble the true orbitals. In contrast the Gaussian 
type-orbitals (GTO) show a wrong shape near and far of the nucleus (no cusp). Far of the nucleus 
the GTOs tend to zero much faster than STOs. 

Nevertheless, the use of STOs was hindered in the last four decades by integration problems 
and the lack of a product theorem as in the case of GTOs. Consequently, Slater orbitals were 
replaced by Gaussian functions in molecular calculations [2] . Despite these difficulties the research 
on Slater orbitals has always continued. In 1982, in a Congress in Tallahassee about Slater orbitals 
Milan Randic described the situation: "Gaussian functions are not the first choice in theoretical 
chemistry. They are used (...) primarily because molecular integrals can be evaluated, not because 
they posses desirable properties. Today this may be a valid reason for their use, but tomorrow 
they may be thought of as bastard surrogates, which served their purpose in the transition period, 
have no longer viable merits and will fail into oblivion" [3]. The need of more accurate molecular 
Configuration Interaction (CI) calculations and Hartree-Fock (HF) calculations of large systems 
has resulted in the use of extremely large basis sets of Gaussian orbitals and dramatically expensive 
calculations. 

The advantage of using Slater orbitals is that a single basis set would consist of one function 
per atom, what would shorten computational times. Also STOs would provide more accurate 
description of molecular properties where the density at the nucleus is important. 

One of the main applications of STOs is the accurate calculation of energy and properties using 
the CI and Hylleraas-Configuration Interaction (Hy-CI) methods [H [5] . For atoms and two-center 
molecules the neccessary integrals are known. It would be desirable to extend the application 
of these methods to three-center molecules. We plan to test the integrals in the case of the H3 
molecule and to use them in CI calculations of the BeH2 molecule. 

As it is well-known there are no analytical general methods to calculate the three- and four- 
center integrals over Slater orbitals. This problem was soon recognized in early years (1959) by 
Mulliken and Roothaan, who called it "The bottle-neck of Quantum Chemistry" [6j. From the 90s 
until today many efforts have been made by several groups to develop efficient algorithms which 
have fructified in new computer programs for molecules. The bottle-neck of these programs is 
the calculation of the three- and four-center repulsion integrals. The approximate methods used 
are: one-center expansion [Tj, [HI El CHI EI], translation [TJl H3], Gaussian expansion [Ti] , Gauss 
transform [T51 [TBJ [T7] , and Fourier transform [TS1 [T91 [201 EB 1221 [23] methods. The disadvantages of 
these methods are on one hand side the low accuracy achieved, what results in longer expansions 
and increased computational times, and on the other side the need of numerical integrations. 
Recently, a special issue of a journal [21] and a book [25] are entirely dedicated to this topic. Also 
improvements of the traslation [13J, Fourier transform [23j . and Gauss transform methods [T7j to 
calculate the three-center repusion integrals have been recently developed. 

The purpose of this paper is to provide the formulas for the analytical evaluation of three- 
center overlap integrals. These integrals are important because the usual two-electron repulsion 
integral of the CI method can be expressed as a linear combination of three-center overlap ones. 
Therefore these overlap integrals should be calculated with high accuracy. 

In 1936 Hirschfelder et al. [2"fi"l |2"T] solved some cases of two-electron repulsion integrals with 
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equal exponents of ls-Slater orbitals in calculations of the H3 molecule. We have extended this 
method to different exponents, which is usually the case in quantum mechanical calculations. 
We will treat here the radial part of the integrals. The angular integration needs to be treated 
separately because the orientation of the angular functions depends on the geometry. Also the 
linear conformation of the centers will be studied first. In Appendix A we give some considerations 
for treating the symmetric triangular case. 

Hybrid-exchange three-center repulsion integrals will be next investigated using this method. 



2. Theory 

We consider a linear three-center molecule, see Figure 1. The distances between the centers may 
be different and R = R a b + R ac . Let b and c be the focii of an ellipse, the atom a is placed along 
R. The elliptical coordinates of electron 1 are defined: 

( r u + ric) (rib - ri c ) 
= ^ » = ^ , 0i = <Pi I 1 ) 

and the volume element is R 3 /8(fi 2 — v 2 )d^,dvdcf). The coordinates of the center a are fi a = 1, 
v a = and <p a = 0. The distance between two points in elliptical coordinates is: 



r\a = ^-(nl + fil + v\ + v 2 a - 2 - - - 1)(1 - u^ifxl - 1)(1 - v 2 a )} 1 ' 2 cos(0i - <j> a )) . 

(2) 

Substituting the coordinates of a in this equation we obtain: 

n a = f(^ + ^ 2 -l) 1/2 , (3) 
The three-electron overlap integral over s-type Slater orbital^] is defined as: 



S na n b n c ((u Cb, Cc R) = j r^e-^r^e-^-'e-^dr, (4) 

where n a ,nb,n c are the principal quantum numbers of the three orbitals located in the centers 
a, b, c, and ( a , (b, Cc are the orbital exponents. In the following we will use the notation r a = r la , 
fb = r ib, t c = ric, Hx = H and V\ — v because the integral depends only on one electron. Expressing 
Tb and r c in elliptical coordinates with focii in b, c, r a is already given by Eq. (3): 

R R 

n = ^(^ + jy), r c =—(n-u) (5) 

Substituting Eqs. (3,5) in Eq. (4), writing explicitly the domains of integration and integrating 
over (p we obtain: 

1 For simplicity only the radial integration will be performed here. 
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x e-^-^g-T^+^-i) 1 / 2 . ( 6 ) 

where the term (/i 2 — z/ 2 ) from the volume element has been already included, a, f3 and 7 are 
orbital exponents whose values are a = R/2(( b + ( c ), (3 = R/2(( b — ( c ) and 7 = R/2( a . Applying 
the Binomial Theorem in the form [28] : 

m+n 

(/i + z/) m (/i - v) n = D™ n (i m+n - k v k , (7) 

k=0 

and defining the coefficients D™ n : 



al(n — a) 
we obtain: 



S nanbnc (CuC b ,C c ,R) = 2na+nb+nc E / ^^e^ 

/ fc=0 71 



^ z/ [Gu 2 + z/ 2 - l) 1 / 2 ] n ^ 1 e-^e-^ 2+l ' 2 - 1 ) 1/2 dz/. (9) 



Now following Hirschfelder, Eyring and Rosen method [26], let us make the change of variable: 

t = (jl 2 + V* - l) 1 ' 2 + U. ( 10 ) 

The limits are transformed as: u = l—^t = /J, + l and for v — — 1— Y t — \i — 1. Defining 
6 = n 2 — 1, we have: 

" = K*-i)' d " = K 1 + ^ dt - (11 < 



Substituting Eq. (11) in Eq. (10): 



(^ 2 + ^ 2 -i) 1/2 = i6 + ?K (12) 



2 V ^ 

and taking into account the transformation of the exponential part of Eq. (9): 

e P" e 711 e vp 7j" e 71 _^ e 2 l 1 J e 7« e 2 e Tig 2 t e 2 e 2 t; (13) 
the three-center overlap integral is: 



4 



o_ T3n a +n b +n c n b+ n c «oo 

w(Ci,c b ,c c ^) = 2WaWrae E / ^-v^ 

! — n J 1 



fc=0 

x / ( 1 + ^) e ^e^-td*, (14) 



t \ t V t 2 



where there is an exponential e 7 * but also one exponential of the type e^ 7 * , difficult to integrate. 
Therefore we will use the following expansion which converges for a few terms of the series: 

fas i A;?! re—, y 

e '•^mT.^{-^ b ) 1 ■ (15) 

iV is the number of terms of the truncated expansion. Eq. (15) is the only infinite expansion 
introduced in the evaluation of the integral. Using the Binomial Theorem: 



< - i) k (* + if 1 + ?) = i (* + ?)"" (* - 9* = S c «"' t "" +M "" (16) 



we get the integral: 



27ri? na+ni,+ric 1 n ^t^ c ^i^J^ TV! nn fc / ' /3 — 7 

Sn a n 6 n c (Cl,C&,Cc,-R) = 2 „ a+rio+nc ^ L L ^fc^" ( ~^T~ 

k=0 u=0 s=0 ' X 

/oo 



x r +1 t n a+k -2u- S -l e -^ dt ^ (17) 



which can be expressed in terms of auxiliary integrals: 

n b +n c n a +k N 



j »6-r«cn a i-re jv 7r_/V! / /3 

fc=0 u=0 s=0 ' ^ 7 

/ (/3 + 7)\ 

X-Bn i ,+n c -fe,s+u;n tl +fe-2u-s-l I - I , (18) 

where B nm .±k(a, (3) are auxiliary integrals, which will be treated in next the Section. 

3. The auxiliary function B nm;± k(ce, f$) 

We have expressed the three-center overlap integrals as a sum of auxiliary integrals B nm .± k (a, (3) 
which are defined: 
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poo /"M+l 

B nm]±k (a,f3)= ^(^-lre-^dfi t ±k e~^dt. (19) 
Ji Jfi-i 

Two cases are distinguished: k is positive t k and the case t~ k . The corresponding undefinite 
integrals are [29] : 

/«. e H» lB = _ e -A^_«_p t >„ (20 ) 

and 

We evaluate now the auxiliary function B nm . +k (a, (3) for > for the actual limits: 

T +1 t k e-^dt = V — ^— /T^ 1 [(/i - l)*- fc e -^-i) _ (/i + ^fe-^-^M+i)] (22) 

Vi h=0 i h - k ) ] 

substituting in Eq. (19), 



5n m;+fc («,/3) = E^Ziyj^" 1 ^/ /i n (/i + l) m (^-l) m+ ^e-(« + ^/i 

- e"' 3 /i n (/i + l) m+fe - fc (// - l) m e- {a+/3)/x rf/i^ . (23) 
The integration over /x can be written as an one-electron auxiliary integral which is well-known: 

k\ 

B nm , +k {a, 13) = 22 (h-k)\ l3 ~ h ~ 1 ( e ^ A n;m,m+k-h{oc + &)- e- p A n . m+k _ Km (a + fij) , (24) 

with 

/oo m+k 
pi n (pi + l) m (pi - l) k e~ a ^ = D™ k A n+m+k ^(a). (25) 
3=0 

In the case of negative k the calculation is messier: 

H+l -fit (_ft\k-l 

dt = ), f \ xt [Et(-(3{fi + 1)) - Ei{-P(jj, - 1))] 



! t k (k-l)\ 

*-l /' ■ ^' r e -/3( M +l) g-^(M-l) 



(ii+l) fc - ft (ii-l) fc - fe 



, (26) 



the auxiliary integral is then: 



B n m;-k(a:, (3) 



-fi) 



k-1 



(k-l)\ 



k-1 



H n {^ 2 - l) m e- a »Ei(-p(n + l))dfi 



/oo 



oo p-/3(M+l) 

//V - l) m e- afl 



k-h d V 



r°° e -/3(M-i) 



. (27) 



simplifying: 



/_/Q\fc-l r roc 

i= o L/i 

/oo 



k-1 

E 



/oo 
/i n (/i + l) m - k+h (ii - l) m e-( a+ ^dn 

/oo 
/i n (/i + !)"> - l) m - k+h e- {a+ ^dfi 



Finally: 



(28) 



(-B) k - 1 m 



i=0 



(-) 



k-1 

-E 



(jfe-/i-l)! 



-Z?)^ 1 [e-*U n;m _ fc+ft , m (a + /3) - e^A n;mim _ fc+?l (a + (3)] (29) 



the auxiliary integral for negative values of A; is a sum of one-electron auxiliary integrals A n;m k(a + 
(3) and a new kind of auxiliary integrals, which will be evaluated in the next Section. 



3.1 Calculation of tL +) («,/3) and T?\a,p) 



The auxilary functions T^} 2m _ 2i (a, (3) and T^ +2m _ 2i (a, (3) need to be calculated. They are defined 
as: 



/oo 
fe-°»Ei(-0(n + l)dii, 

/oo 
li n e-*»Ei(-P(ii-i)dii. 



(30) 



The auxiliary functions Tn + \a,f3) can be calculated using recursion relations. The first term is: 



/oo 
e- a »Ei(-p(n + l)dn. 

The higher terms can be calculated using the differential equation: 



Ti +) («,/3) = (-ir^ (+) (a,/3). 



Let us first solve T^ + \a, /3), making the change of variables: 



y 



P(n + l)=y, /3dfi = dy, n=--l, 

the integration limits of the domains are transformed u = 1 — > y =4> 2(3 and u — > oo 
To calculate Tq + \cx,/3) one has to solve the integral: 



-e<y i )Ei{-y)dy. 



roc -i 

T (+ W) = / 

J2/3 P 

Using the integration formula: 

/OO 1 1 
-e-^Ei(-y)dy = - [Ei(-c)e^ c - c£i(-( 7 + 1))] , 



we have 



T (+ W) 



e a 1 



Ei(-2.i)r-p p - Ei[ -[^ + 1)21 



a 



and simplifying: 



The nth-term is: 



T (+) (a, 0) = —Ei(-2p) - —Ei(-2(a + /?)). 
a a 



d n / p~ a p a 
Ti +) (a,(3) = (-1)"^ l—Ei(-2p) - -Ei(-2( a + fi)) 



Let us use the Leibnitz formula: 

d n 



dx r 



" Jn—m Am 

(m g ( X )) = ^c m ^f(x)— g (x) 



m=0 



'dx r 



dx % 



calculating £,^Ei(-2p) and Ei(-2(a + /?)): 



Jn „-o; 

{ - l)n d^^ El{ ~ m = ( " 1)n ^ (^o(«)^(-2/3)) = M«)Ei(-m 



da T 



(31) 
(32) 

(33) 
y — > oo. 

(34) 
(35) 



(36) 



(37) 



(38) 



(39) 



(40) 



and 
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Jn p a Jn 

{ ~ ir d^a Et{ ~ 2{a + = { ~ ir d^ (-M-a)Ei{-2{a + /?))) 

= (-^J2C^M- a) —E t{ -2( + m 

m=0 

n 

= -A n (-a)Ei(-2(a + /?)) + ^(-l)"- m ^A n _ m (-a)^(-2(a + /?)) 

m=l 

We have used the following definitions of the exponential integral: 

pz e -t poo e -t 

Ei(-z) = / — dt = - / — dt 

1 J J-oo t J z t 

f°° e _t 

E x {-z) = / — dt = -Ei(-z) = Ei(l, z). 

J z t 

The derivatives of the exponential integral are 



/oo 
e -^+W dt = 2A (2(a + f3)), 



d 

da 

d m , rf m ' ! 



£i(-2(a + fi)) = — =T 2A (2(a + /?)) = (-l)— 1 2 m A m _ 1 (2(a + /?)) 



(io; m da 

Jn f a \ 

(-1)"^ ^^(-2(a + /3))J = -(-irA t (-«)^(-2(a + /3)) 

+ ^(-l)«-™2 m C"A n _ m (-a)(-l) m - 1 A m _ 1 (2(a + /?)) 

m=l 

Finally we have: 

Ti + \a,P) = A n (-a)Ei(-2f3) + (-l) n A n (-a)Ei(-2(a + /3)) 

n 

- ^(-l)"- m 2 m C™ A n _ m (-a)(-l) m - 1 A m _ 1 (2(a + /?)) 

m=l 

To calculate the auxiliary functions T^~\a,0) defined as: 

/oo 
fM n e-^Ei(-0(jM - l)dn 

we use also recursion relations. Let us calculate also Tjj \oi,0) defined as: 

/oo 
e- a "Ei(-0(» - l)dfi. 



The following change of variable is done: 



u = Ei(-P(n-1)), ^du= ^— — 



dv = e" a/ \ 

a 

The integral Tff~\a, /3) is the limit: 



e 



/■oc 

T {_) (a,/3) = lim / e~ a ^(-/3(/i - l)d/x, 

integrating by parts: 

Tt\a,P) = lim (—e-°*Ei(-l3(iM-l) + - f 



i r°° p -/3(m-i) 



making again a change of variables, in this case: /i — 1 = u ^> dfi = du and /x = 7 =>• m 
with /i — )> 00 =>- u — > 00, we have: 

/ 1 p-a /-oo -( a +/3)u \ 

T ( ~>, /3) = lim e - Q ^(-/3(/i - 1) + / du 



e 



'7-1 

-a 



lim (Ei(-P(fi - 1) - £i(-(a + - 1)) 



00 e -t ,oc e - t 



a 7->i 

lim ( - / — di f / : — f// 

using the identity 



7^1 V Jfa-l) t i (Q+/3 )( 7 _i) t 



r e -dt=r{ i -^\dt-r e -dt-C 

Ja t J a \ t J h t J a 



dt 
T 



( ' = I I — - — J dt- J —dt = 0.5772156649015. 



and the definition of the Euler constant C: 



we get: 



T ( - } (a, 0) = — lim <| / 1 ^^rft - f° — di - f 1 - - f°° l -^—dt 



a [J 13(^-1) t Jl t i/3(7-l) * J(a+p)(y-l) * 



/ rft+ / — 

il * 7(a+/3)(7-l) * 



10 



Tt\a,(3) 



a 7->i 



lim [C-C- hit 



1 

/3(7-l) 



hit 



1 

(a+/3)(7-l) 



(57) 



simplifying: 



T ( - } (a,/3) 



— lim (ln[/3( 7 - 1)] - ln[(a + /3)( 7 - 1)]) 



— lim (In (3 + ln( 7 - 1) - In (a + (3) - ln( 7 - 1)) (58) 
a 7->! 



we finally get for Tq ^(a, /3) a simple expression: 



T Q ( ~W) 



In 



P 



a \a + (3 



Now we consider the nth-term: 



(-ir-^T «(a,/3) = (-1)^— In ( -A- 



-r 



A (a) In 



/3 



a + (3 



n Jn—m Jm / o 



m=0 



a + (3 



obtaining finally: 



Ti-\aJ)=A n (a)\n 



where A n (a) is the Mulliken integral: 



P 



a 



P 



m=l 



im 



(a + py 



-A. 



la 



(59) 



(60) 



(61) 



AJat) 



t n e~ at 



me sr-^ 

a n+l 



fc=0 



a 



(62) 



In Table 1, the pattern of convergence of the three-center overlap integral is shown. To obtain 
a convergence of 15 decimal digits 19 terms of the expansion are needed. To obtain 30 decimal 
digits accuracy, about 30 terms of the expansion are required. These values have been calculated 
with Maple [30]. This algorithm will be programmed in a Fortran source code and the efficiency 
of the calculation of the integral can be tested. Further, this integral will be tested by comparing 
the two-electron three-center integrals with the values obtained by the SMILES computer code 
[31] version which uses the Gaussian expansion method. 

In Table 2, a number of three-center one-electron integrals for different exponents, quantum 
numbers and interatomic distances is presented, with respect to the number of terms neccessary 
in the expansion. 
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3.2 Summary 

The three-electron overlap pone-electron integral is: 



Sn a n h n c (Cl? (b, Co R) 



n b +n c n a +k N / , 

_ V V V 7FiV! D n b n cD n a k ( P ~ 7 

fc=0 m=0 s=0 



X Bn b +n c — k,s+u;n a +k— 2u— s— 1 I 



The auxiliary integrals for postive and negative k: 



(63) 



fc! 



5 nm;+fc (tt, 0) = (frjfc)^ ^ 1 (^A.;m,m+fc-/i(a + 0) - e ^A.;m+fc-h,m(Q! + /?)) (64) 



/i=0 



i=0 



^n+2m-2i( a i P) ^n+2m,-2i( a ^ P) 



fe-1 

-E 



(k-h- 1)! 



The auxiliary functions: 



-/J)* -1 [e- p A n . m _ k+h)m {a + p) - e^A n „. k+h {a + (3)] (65) 



T«(a,/3) = A n (-a)Ei(-2f3) + (-l) n A n (-a)Ei(-2(a + f3)) 

n 

_ ^(_ 1 )-2T;i_ m (- a )(-l) m - 1 A m _ 1 (2( (l + /?)) 



m=l 



^(a,)?) = -Un) in ( ; + £ r;;. ( "' ,„(,, j 



m=l v ' 



The one-electron auxiliary integrals: 

A n (a) = 



t n e~ at 



m+k 



n\e a v-^ a k 



a 



fcf' 



n+1 fc| 

fc=0 



Ai;mfc(aO = ^ Df k A n+m+k -j (a) . 
3=0 



(66) 



(67) 



(68) 
(69) 
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4. Two-electron three-center integrals 



There are several kinds of two-electron integrals where three centers are involved, see Figure 1: 
Coulomb-exchange (the first two) and hybrid-exchange (the last two): 

[aa,bc} = J J Xa(l)Xa(l)^X6(2)x c (2)rfr 1 rfr 2 , 
[bb,ac] = J J X6(l)X6(l)^Xa(2)x c (2)dr 1 dr 2 , 
[ab,ac] = J J Xa(l)X6(l)^Xa(2)Xc(2)rfridr 2 , 

[ab,bc] = j J Xa(l)x b (l)^-Xb(2)Xc(2)dr l dr 2 . (70) 

The Coulomb exchange integrals can be evaluated by integration over the coordinates of electron 
1 (32], obtaining a linear combination of three-center overlap integrals calculated in this paper: 
In general a two-electron integral can be written: 



oo roo 



[aa,bc]= dr 2 Xa(l)Xa(l)gi2Xb(2)Xc(2)dT 1 dT 2 , (71) 

JO Jr 2 

where g\ 2 is an operator of r\ and r 2 . For L = the operator is g\ 2 = g\ 2 where g\ 2 means the 
largest of r± and r 2 . There are two integration domains: r% < r 2 and T\ > r 2 . One may refer the 
interelectronic distance r± 2 to center a and expand it in terms of the variables r\ a and r 2o : 



roc roc j roo roo j 

[aa,bc}= / dn / Xa(X)Xa(X)— Xb(2)Xc(2)dr 2a + / dr 2a / Xa(l)x„(l) — Xb(2)x c (2)dr la 

JO Jri a r 2 JO Jr 2a T ^ 

koo rr 2a r<x> r<x> i 

dr 2a X a(l)Xa(l)—Xb(2)Xc(2)dr la + dr 2a X a(l)Xa(l)—Xb(2)Xc(2)dr la . (72) 

Jo r 2a Jo Jr 2a r la 

As these orbitals are located at the same center, we can integrate over the coordinates of one- 
electron [32J: 

A n (x, a) = j r? a e-^dr la = e — £ (73) 



X 



k=0 



/x | | n k 

rte—dr la = A n (a) - A n (x, a) = — - e — £ T^MT^m ' ( 74 ) 

The Coulomb-exchange integrals of Eq. (70) can be evaluated integrating over the coordinates of 
electron 1, writting them explicitly: n\ a — n\ +rii — 1 and C = Ci + Ci : 

[aa,bc]=A7i J r2~ 1 e-^ b r^- x e-^ nc dT 2 J r^V^ 1 " ^-r 2 la dr la , 

[66, ac] =4tt ^ r^e-^r^e'^d^ J r^e'^ ^r 2 lb dr lb , (75) 
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using Eqs. (73) and (74) we obtain 



poo f i pria poo ~\ 

[aa,bc]=A7r V^V^rV^dra \ — / rZ +1 e~ Caria dr la + r^e~^ na dr la \, 

J0 l r 2a JO Jr 2a ) 

(76) 

[aaM - - r*w^ { - 



T2a Ca na+2 ^2a C a " 

_ ^ (n a + l)\ n a -k -( a r 2a , ( n «) ! „n a -fc -Car 2a I / 77 ^ 

^ K + l-A;)!^^ + ^K-A;)!e ir2a f' U ° 

aa, 6c] = 4vr /°° r^V^V^V^Ta ( — (l - e^ r2 °) 

fc-1 K)! 
K + l-A:) (n a -A;)!C* 



+ V - - ^ r ^-k -QaV 2a I / ?8 ^ 



k=0 

Finally the Coulomb-exchange reulsion integral is expressed as a linear combination of three- 
center overlap integrals defined in Eq. (4): 



f na _|_ i)l 

[aa, be] = Ait — - ' (S Qjnb>nc (0, ( b , ( c , R) - S , nb ,n c {( a , Cb, Cc R)) 



+ 4VT g K + ! _ fc) K _ a k)lQ+1 Sn M n b ,n c (Ca, Cb, Cc, A), (79) 

where the integral So t n btric (0, Cb, Cc -R) is the three-center nuclear attraction integral defined in Ref. 
[33]: 

dr (80) 

fa 

So,n b ,n c ((a, Cb, Cc, R)) should be analyzed. Anagously for [bb,ac] we obtain: 



[66, ac] = 47T ^ b n "^ + ^' (5 ,„ a ,„ c (0, Co, Cc R) ~ S ,n a ,n c {(b, Ca, Cc R)) 

Cb 



"6 



+ 47r E („ 6 + i_jfc) K-fc)!C b fc+1 ^- fc+1 ' na ' nc(Cf " Ca ' Cc ' ^ (81) 

,n. a ,n c (0, Ca, Cc, -R) are expressed as: 

dr (82) 

r b 

The integral So )na>nc (Cb, Ca, Cc, -R) is an especial case of overlap integral with n a = —1 and should 
be treated separately, as Eqs. (9-14) are derived for n a > 1. 

Next we will calculate some two-electron integrals and compare with the literature values 
and values obtained with a computer program by Rico et al. [3T] using the Gaussian expansion 
method. 
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Conclusions 



In this paper three-center overlap integrals are calculated, which are used to calculate the two- 
center three-center repulsion integrals. The integrals consist of an expansion which converges 
relatively fast. 
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Appendix A 

In Figure 2 the symmetric triangular conformation of the centers is shown. The interelectronic 
distances are equal R = R ab = R ac . The elliptical coordinates of the center a with respect to the 
focii b, c are fixed values: 

Rab ~\~ Rac n / A 1 \ 

Ha = ^ = 2, v a = 0, (A.l) 

and the relation: 

cos(0i — <f> a ) = cos 0i cos a + sin 0i sin a (A. 2) 

independently of the choose of the angle a , in integrations over s-type Slater the integrals 
J Q 27r cos0i(i0i = and J Q 27r sin 0i<i0i = vanish. Then we have: 

r\ a = + + (A.3) 

and therefore 

r la =^\ + ul + 2fl\ (A.4) 

accordingly the substitution in this case is: 

t = (/i 2 + v 2 + 2) 1 / 2 + v. (A.5) 

which is similar to the one used in Section 1, and therefore the steps of integration described in 
previous Sections can be used. The triangular conformation is very important because it is the 
ground state conformation of the H 3 molecule and many molecular systems. 
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Figure 1: Definition of the coordinates of one electron with respect to three centers: linear 
conformation. 

Figure 2: Definition of the coordinates of one electron with respect to three centers: triangular 
conformation. 
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Table 1: Convergence pattern of the three-center overlap integral for the orbital parameters 
n a — 1, rib — 2, n c — 1 and ( a = 1.6, 0> = 1-4, Cc = 1-2; -R = 1.4 a.u. 



1 


-0. 


11416 










2 


-0. 


12293 


3 








3 


-0. 


12160 


22 








4 


-0. 


12219 


350 








5 


-0. 


12207 


7743 








6 


-0. 


12210 


04107 








7 


-0. 


12209 


67757 


3 






8 


-0. 


12209 


73204 


83 






9 


-0. 


.12209 


72448 


624 






10 


-0. 


.12209 


72547 


7065 






11 


-0. 


.12209 


72535 


25580 






12 


-0. 


,12209 


72536 


75642 







13 


-0. 


12209 


72536 


58049 


61 




14 


-0. 


12209 


72536 


60052 


409 




15 


-0. 


12209 


72536 


59828 


9497 




16 


-0. 


12209 


72536 


59853 


35836 




17 


-0. 


12209 


72536 


59850 


73400 


8 


18 


-0. 


12209 


72536 


59851 


01158 


63 


19 


-0. 


12209 


72536 


59850 


98261 


802 


20 


-0. 


12209 


72536 


59850 


98559 


9711 


21 


-0. 


,12209 


72536 


59850 


98529 


65861 


22 


-0. 


,12209 


72536 


59850 


98532 


70074 8 


23 


-0. 


,12209 


72536 


59850 


98532 


39925 80 


24 


-0. 


12209 


72536 


59850 


98532 


42873 615 


25 


-0. 


12209 


72536 


59850 


98532 


42589 4768 


26 


-0. 


12209 


72536 


59850 


98532 


42616 42719 


27 


-0. 


12209 


72536 


59850 


98532 


42613 91777 


28 


-0. 


12209 


72536 


59850 


98532 


42614 14624 


29 


-0. 


12209 


72536 


59850 


98532 


42614 12602 


30 


-0. 


.12209 


72536 


59850 


98532 


42614 12774 



Table 2: Three-center overlap integrals. N is the number of terms of the expansion needed to 
calculate accurately the printed number of decimal digits. 



n a 


n h 


n c 


(a 


Cb 




R 


Sn a n b n c {C,li 0>> Cc? R) 


N 


1 


1 


1 


1.3 


1.3 


1.3 


1.0 


-0.10852 96351 53609 38988 


25 


1 


1 


1 


1.1 


1.3 


1.5 


1.0 


-0.41057 77481 70340 63986 


26 


1 


2 


2 


2.0 


2.0 


2.0 


2.0 


0.00409 87551 44074 883 


30 


2 


2 


2 


1.6 


1.4 


1.4 


2.5 


0.08943 10170 26917 0505 


28 


4 


3 


2 


2.6 


2.4 


1.6 


3.0 


0.00362 95811 14392 3873 


25 



19 



